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The adiabatic selfconsistent collective coordinate method is applied to an exactly solvable 
multi-C'(4) model which simulates nuclear shape coexistence phenomena. Collective mass 
and dynamics of large amplitude collective motions in this model system are analysed, and 
it is shown that the method can well describe the tunneling motions through the barrier 
between the prolate and oblate local minima in the collective potential. Emergence of the 
doublet pattern is well reproduced. 



§1. Introduction 



in 
o 
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\ Microscopic description of large amplitude collective motion in nuclei is a long- 

standing fundamental subject of nuclear structure physics.^''^ In spite of the steady 
development in various theoretical concepts and mathematical formulations of them, 
CJ ' application of the microscopic many-body theory to actual nuclear phenomena still 

2 ■ remains as a challenging subject^^"^''^ (see Ref. 1311) for a recent comprehensive re- 

view) . Shape coexistence phenomena are typical examples of the nuclear large ampli- 
tude collective motion and have been investigated from various points of view'^'^^"^^) 
^ ■ For instance, even in typical spherical nuclei like Pb and Sn isotopes, excited de- 

, formed states have been systematically observed in low-energy regions,'^ and the 

coexistence of prolate, spherical and oblate shapes have been recently reported for 
i86,i88p]2_EHl As another example, we mention a recent discovery of coexisting two 
rotational bands in ^^Se, which are associated with oblate and prolate intrinsic 
shapes.!^ These are only a few examples among abundant experimental data. 
Clearly, these data strongly calls for further development of the theory which is able 
to describe them and renews our concepts of nuclear structure. From the viewpoint 
of the microscopic mean-field theory, these phenomena imply that different solutions 
of the Hartree-Fock-Bogoliubov (HFB) equations (local minima in the deformation 
energy surface) appear in the same energy region and the nucleus undergoes a large 
amplitude collective motions connecting these different equilibrium points. Identi- 
ties or mixings of these different shapes are determined by the dynamics of such 
collective motions. 

On the basis of the time-dependent Hartree-Fock (TDHF) theory, the selfcon- 
sistent collective coordinate (SCC) method was proposed as a microscopic the- 
ory of such large amplitude collective motions.^' It was extended to the case 
of time-dependent HFB (TDHFB) including the pairing correlations and has 
been successfully applied to various kinds of anharmonic vibration and high-spin 
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rotations^^^"^''^ In order to apply this method to shape coexistence phenomena, 
however, we need to further develop the theory, since the known method of solving 
the basic equations of the SCC method, called ?7-expansion method,!^ assumes a 
single local minima whereas several local minima of the potential energy surface 
compete in these phenomena. Quite recently, a new method of solving the basic 
equations of the SCC method, called adiabatic SCC (ASCC) method, has been pro- 
posed P This new method uses an expansion in terms of the collective momentum 
and do not assume a single local minima, so that it is expected to be suitable for a 
description of the shape coexistence phenomena. The ASCC method may also be 
regarded as a successor of the adiabatic TDHF (ATDHF) methods. It inherits the 
major achievements of the ATDHF theory (as reviewed in Ref.l31|l) and, in addition, 
enables us to include the pairing correlations selfconsistently removing the spurious 
number fluctuation modes. 

The major purpose of this paper is to examine the feasibility of the ASCC 
method for applications to actual nuclear phenomena. This is done by applying 
it to an exactly solvable model called multi-0(4) model and testing results of the 
ASCC method against exact solutions obtained by diagonalizating the Hamiltonian 
in a huge bases. This solvable model may be regarded as a simplified version of the 
well-known pairing-plus-quadrupole (P+Q) interaction model.'^''^ Namely, only 
the K = component of the quadrupole deformation is considered in a schematic 
manner, and the model has been widely used as a testing ground for various mi- 
croscopic theories of nuclear collective motion.'^''^''^J^ The multi-0(4) model 
possesses a symmetry, analogous to the ordinary parity quantum number, with re- 
spect to the sign change of the "quadrupole" deformation. Accordingly, it can be 
utilized as a simple model of many-body systems possessing double well structure 
where large amplitude tunneling motions take place through the barrier between the 
two degenerate local minima of the potential (which correspond to the prolate and 
oblate shapes). Because of the special symmetry of the model, the "prolate" and 
"oblate" shapes mix completely. Of course, in contrast to the ordinary parity, such 
an exact prolate-oblate symmetry does not occur in reality, and in this sense this 
solvable model is somewhat unrealistic. Nevertheless, using this model, we can make 
a severe test of the theory by examining its ability to describe the "parity doublet" 
pattern. 

In this paper, we would like to focus our attention on the collective dynamics 
and the collective mass of large amplitude collective motions. Needless to say, the 
barrier penetration depends on the collective mass in a quite sensitive manner. A 
similar investigation on the collective dynamics of this model was done in Ref . I27|) , 
in which the number fluctuation's degree of freedom was explicitly removed from the 
model space. In the present approach, the spurious number fluctuation modes are 
automatically decoupled from the physical modes within a selfconsistent framework 
of the TDHFB theory. This will be a great advantage when the method is applied 
to realistic nuclear problems. 

This paper is arranged as follows: In § 2, the basic equations of the ASCC 
method are recapitulated. In § 3, a brief account of the multi-0(4) model is given. 
In § 4, we apply the ASCC method to the multi-0(4) model and derive explicit ex- 
pressions necessary for numerical calculations. In § 5, we present results of numerical 
analysis, and conclusions are given in § 6. 
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§2. Basic equations of the ASCC method 

In this section, instead of recapitulating the general outline of the SCC method, 
we summarize the adiabatic version of them formulated in Ref. I58|l . We assume 
that large-amplitude collective motions are described by a set of the TDHFB state 
vectors \4>{q,p,(p, N)) that are parametrized by a single collective coordinate q, the 
collective momentum p conjugate to q, the particle number and the gauge angle 
if conjugate to A''. The time evolution of \(j){q,p,ip, N)) is determined by the time- 
dependent variational principle 

6{4>{q,p,ip,N)\i-^^-H\4>iq,p,^,N)) =0. (2-1) 
As discussed in Ref. , we can set 

\^{q,p,ip,N))=e-'''^\cP{q,p,N)). (2-2) 

Then, since the Hamiltonian H commutes with the number operator A^, the gauge 
angle (p becomes cyclic. The basic equation of the SCC method consists of the in- 
variant principle of the TDHFB equations ()2-l() and the canonical variable condition 
which can be written as equations for the state vectors \<j){q,p, N)) , 



{<P{q,p,N)\i—\cl){q,p,N))=p, (2-3a) 

{4>{q,p,N)\^\(t>{q,p,N)) =0, (2-3b) 
top 

{(l){q,p,N)\N\(P{q,p,N)) =N, (2-3c) 

{HQ,P,N)\4K^\<P{q,p,N)) =0. (2-3d) 



The third equation guarantees that the particle-number expectation value is kept 
constant during the large-amplitude collective motion described by the collective 
variables {q,p)- Assuming that the large-amplitude collective motion described by 
the collective variables {q,p) is slow, we now introduce the adiabatic approximation 
to the SCC method. Namely we expand the basic equations with respect to the 
collective momentum p. Since the particle-number variable is a momentum vari- 
able in the present formulation, we also expand the basic equations with respect to 
n = N — No, when we consider a system with particle number A'^o- We then keep 
only the lowest order term. The TDHFB state vectors are thus written as 



|,/.((?,p,Ar)) = e*P«(«)+-^(«)|0(g)), (2-4) 

where Q{q) and 0{q) are infinitesimal generators with respect to {(piq))- We also 
define an infinitesimal generator P{q) by 

e-W<?)|^(^)) = |^(g + 5g)). (2.5) 



We insert the TDHFB state vectors ()2-4() into ()2-l() and make an expansion with 
respect to p and n. Requiring that the time-dependent variational principle (the 
canonical variable condition) be fulfilled up to the second (first) order, we obtain 
the basic set of equations of the ASCC method to determine the infinitesimal gen- 
erators Q{q) and P{q) as follows: 
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Canonical variable conditions 

{Hq)\[QiQ),HQ)]\HQ))=i, (2-6a) 
{cPiq)\mQ),mmq))=i, (2.6b) 

the other expectation values of commutators among {Q{q),P{q),0{q),N} being 
zero. 

HFB equation in the moving frame 

6{cj,{q)\HMiq)\HQ))=0, (2-7) 
where the Hm{q) is the Hamiltonian in the moving frame defined by 

HM{q)=H-Xiq)N-^Q{q). (2-8) 

Local harmonic equations 

6 [HMiq], Qiq)] - \B{q)P{q) = 0, (2-9a) 

5 {<i>{q)\ [HM{q), ]p{q)]-C{q)Q{q)-^^[[HM{q), {H-X{q)N)A],Q{q)]-^N \<P{q)) = 0, 

(2-9b) 

where the local stiffness C{q) is defined by 

^/ , d^V 1 dBdV , , 

and {H — XN)a indicates the a^a^ and oo parts of {H — XN) containing two- 
quasiparticle creation and annihilation operators. The collective potential V{q), 
the inverse mass B{q), and the chemical potential X{q) are defined below. Equa- 
tions 1)2 -yal) and ()2-9b|) are linear equations with respect to the one-body operators 
Q{q) and P{q)- They have essentially the same structure as the standard RPA 
equations except for the last two terms in Eq. ()2-9b|) , which arise from the curvature 
term (derivative of the generator) and the particle-number constraint, respectively. 
The infinitesimal generators Q{q) and P{q) are thus closely related to the harmonic 
normal modes locally defined for \(j){q)) and the moving frame Hamiltonian HMiq)- 
The collective subspace defined by these equations will be uniquely determined once 
a suitable boundary condition is specified. 
The collective Hamiltonian is given by 

n{q,p, N) ^ {(t>{q,p, N)\ H\(t>{q,P, N)) (2-lla) 

= V{q) + ^B{q)p^ + X{q)n, (2-llb) 

up to the second order in p and the first order in n, where 

V{q) = n{q,p, N)\,=o,N=No = mQ)\ H \Hq)) , (2-12) 
i a^H(g,p ,iV) 
2 5p2 



Big) = ^ \p=o,N=No = - {<Piq)\ [[H, Q{q)],Qiq)] \Hq)) , (2-13) 



M'/) = ^^^^^'^^ =0.^=^0 = {4>iq)\ [HMq)] mq)) ■ (2-14) 
For the system with N = Nq particles, we can put n = 0. 
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§3. Multi-0(4) model 

The multi-0(4) model may be regarded as a simplified version of the con- 
ventional P+Q interaction modelj^''^ where only the K = component of the 
quadrupole deformation is considered in a schematic manner. It has been used for 
schematic analysis of anharmonic vibrations in transitional nuclei and of various 
kinds of large-amplitude collective motion .'^^'^''^''^ 

We define bilinear fermion operators for each j-shell, 

= X/ 4'n4-m' Bj ^ X/ (3-la) 

m>0 m>0 

Nj = cl^Cjm, Dj = "^jrnC^^Cjm, (3-lb) 

m m 

with 

_r 1 \m\<f2,/2 

~\ -I \m\ > [2j/2. ^'^ ^' 

The sign of ajm is chosen so as to simulate the behavior of the quadrupole matrix 
elements {jm\ r^l20 |j"^)> and we assume that the pair multiplicity I7j = j + ^ is 
an even integer. The set of operators {^], Aj, Bj, Bj, Nj, Dj} form a Lie algebra of 
0(4). We then define their extensions to the multi j-shell case, 

^t = ^At, B^ = Yb], N = Yn„ b = Yd,D,, (3-3) 

3 j j j 

the coefficients dj in D simulating the magnitudes of the reduced quadrupole matrix 
elements of the j-shells, and introduce a model Hamiltonian 

H = ho- ^GiAU + AA^) - ]^xb\ (3-4a) 

i 

where denote single-particle energies of the j-shells, and G and x represent the 
strengths of the pairing and the "quadrupole" interactions, respectively. Note that 
this Hamiltonian is invariant with respect to the conversion of single-particle states 
(j, ±|m|) ^ {j,±{Qj — \m\)). Thus, eigenstates can be classified according to the 
"parity" quantum number associated with this symmetry. 
If we make different combinations of these operators as 



Kj+ = 


^(^] + ^])' 


Lj+ - 






(3-5a) 


Kj- = 


^iAj+Bj), 


h- = 




-Bj), 


(3-5b) 


KjO = 




Ljo = 






(3-5c) 



the sets {Kj+, Kj^, Kjo) and {Lj+, Lj^, Ljo) separately form SU{2) algebras, and 
they commute with each other. Namely, the multi-0(4) model is equivalent to the 
multi SU{2) SU{2) model. Thus, we can diagonalize the Hamiltonian (|H-4j) in a 
basis set 

J\\nK,,nL,) = X{{K,+r<^{Lj+r-m. (3-6) 
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to get the exact eigenvalues and eigen-vectors, where nxj and nij respectively indi- 
cate numbers of the K and L pairs in the j-shell. They satisfy < nxj, nLj < ^jf^ 
and ^j{nKj + n^j) = No/2. We note that, in the special case that single-particle 
levels are equidistant, all dj are equal, and all f2j = 2, this model reduces to the 
one used in Ref . ^5) which studied collective mass in finite superconducting systems. 

§4. Application of the ASCC method to the multi-0(4) model 

4.1. Quasiparticle representation 

We are now in a position to apply the ASCC method to the multi-0(4) model. 
For separable residual interactions such as those in this model, it is customary to 
neglect their Fock terms.'^''^ We follow this prescription. Accordingly, in the 
following, we use the notation HB in place of HFB. It is readily seen that the TDHB 
state vectors \(p{q)) in the multi-0(4) model can be written in the BCS form, 

\<l>{q))=exp{Y^e,{q)iAl-A,)}\0), (4-1) 

i 

where 9i{q) are related with the coefficients Ui and Vi of the Bogoliubov transforma- 
tion to the quasiparticle operators a| and Oj, 

as Ui = sin 6i and Vi = cos^j. Here, i = {j,m), —i = {j, —m), and denotes a sum 
over levels with m > 0. We use these conventions hereafter. 

The pair operator = clc^_j^ and the number operator Ni = (^^Ci+c^_^C-i for the 
degenerate single-particle levels (i, —i) are written in terms of the quasiparticle-paiv , 
A| = a|alj and quasiparticle-numhei operator, iVj = a|aj + a^ja_j, as 

a] = UiVi + u^aI - vfAi - UiViNi, (4-3a) 
= 2vf + 2uiVi{Al + Ai) + (nf - vf)Ni. (4-3b) 

The quasiparticle operators, A|, Aj, and iVj, satisfy the commutation relations 

6^A1-N^), (4-4a) 
26u>aI. (4-4b) 

We define the deformation D{q) and the pairing gap A{q) for the TDHB state {(piq)) 
as 

Dig) = ((^(g)|D|<^(g)) =2Y,diaivl (4-5) 

i 

A{q) = At \cPiq)) = G^um. (4-6) 

i 

4.2. Quasiparticle RPA at local minima 

We start from the standard procedure for describing small-amplitude vibrations 
around the local minima of the collective potential. Namely, we apply the quasipar- 
ticle RPA about the HB equilibrium points. Since the ASCC is equivalent to the 



A- 4"^ 



Ni,Al 
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HB+RPA at equilibrium states, the quasiparticle RPA modes provide the boundary 
condition for solving the local harmonic equations of the ASCC method. 

For the equilibrium HB state |i;^o)) the HE equation is given, as usual, by 

^"o'° -.f: .„)(:: 

where Aq and Aq denote the pairing gap and the chemical potential, and 

ei = e- - xdiCTiDo (4-8) 

are single-particle energies at the equilibrium deformation Dq. The quasiparticle 
energy and the particle number are written as Ei = \J (cj — Aq)^ + = 

Writing the RPA normal coordinates and momenta as 

i 

Pn = iY^Pl[A\-A,), (4-9b) 

% 

we can easily solve the quasiparticle RPA equations, 

b [H - XoN, Qn] - -BnPn \M = 0, (4-lOa) 

I 

6 {^o\ [H - XoN, -Pn] - CnQn ^ = 0, (4- 10b) 

I 

where Bn and On denote the inverse mass and the stiffness for the n-th RPA so- 
lution. Note that the local harmonic equations (|2-9aj) and ()2-9b|) reduce to the 
RPA equations at the HB equilibrium points, since the third and fourth terms in 
Eq. ()2-9b|) vanish there. The RPA dispersion equation determining the frequencies 
w„ = VBnCn is given by 

detiSkk') = (4-11) 
where S" is a 3 x 3 matrix composed of 



Sn 


= Y^iff/m - 

i 


■uj)+gl/{2Ei+u:))-l, 


(4-12a) 


Sl2 


= Y,{fi9i/m 

i 


-u:)+gJi/{2Ei+io)), 


(4-12b) 


Sl3 


= J2{fihi/{2Ei 

i 


-u;)+gihi/{2Ei+uj)), 


(4-12c) 


S22 


= Y.{gVm - 

i 


u;)+ff/{2E,+uj))-l, 


(4-12d) 


S23 


i 


-uj) + nhi/{2E,+uj)), 


(4-12e) 


S33 


= Y.{hy{2E, - 


uj) + h^/{2E,+uj))-l, 


(4-12f) 



and Sk'k = Skk', with the notations fi = VCuf, gi = VCvf, and hi = 2^/xdiaiUiVi. 
If X = 0, the above dispersion equation reduces to 
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which involves two weh- known quasiparticle RPA normal modes; the pairing vibra- 
tion (w ~ 2 A) and the pairing rotation (uj = 0). On the other hand, if we neglect 
the residual pairing interactions, it reduces to 

which involves a normal mode analogous to the (3 vibrations in deformed nuclei. We 
note that these three kinds of normal modes are decoupled at the spherical point 

{D = 0) where (0(g)| [z), + a\ \^{q)) = 0. 

4.3. HB and local harmonic equations in the moving frame 

In order to find a collective subspace in the TDHB space, we need to solve 
the RPA-like equations in the moving frame; the local harmonic equations. These 
equations determine generators of the collective space, Q{q) and P{q). Now let us 
summarize the local harmonic equations for the multi-0(4) model. With notations 
F,=i = A, Fs=2 ^ -D, and 

Fi±)^(F,±F])/2 = ±Fi±)t, (4-15) 
the multi-0(4) Hamiltonian is written as 

H = h,-\Y: ^.fWfW + \ E -sF^~'^n-\ (4-16) 

s=l,2 s=l,2 

where the suffices s = 1 and 2 indicate the pairing and the "quadrupole" parts, 
respectively, and ki = 2G and K2 = X- The equation of motion for the time- 
dependent mean-field state vector |0((j',p)) is written as 

H<l^{q,p)\i^^-hm,p)) = ^, (4-17) 
with the selfconsistent mean-field Hamiltonian 

h = ho-Y, ^sH-^^ {<P{q,p)\ m,p)) + J^k.F]-) {<^{q,p)\ Fi") \<^{q,p)) . 

s s 

(4-18) 

The HB equation in the moving frame (|2-7|) and the local harmonic equations (|2-9ajl - 
1)2 -Obi) then become 

5m)\hM{q)\<P{q)) = ^, (4-19a) 
Hm\ [hM{q),Qiq)]-Y,fi^'!F(-^ - ^B{q)Piq) = 0, (4-1%) 

s 

S {Hq)\ [hdq) , ]B{q)P{q)] - J2 fi^^F^ _ B{q)C{q)Q{q) - J] f^+}F^+^ 

s s 

- Y.[F^-\ iHq) - X{q)N)A]f^Q-! - JnN = 0, (4.19c) 
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where /im(9) is the selfconsistent mean-field Hamiltonian in the moving frame de- 
fined by 

dV 

hM{q) = Hq)-X{q)N -—Q{q), (4-20a) 
Kq) = ho-Yl '^^^i^^ ^s"-^ l<^(9)) ' (4-20b) 



and 



fq,! = -^s m)\ [F^-\Q{q)] \m) , (4-21a) 
/JJ) = {m\ \B{q)P{q)] \m) , (4-21b) 



fR.! = -n'^s im [[Fr\CKq) - mN)A],Q{q)] , (4-21c) 



f^ = B{q) — . (4-21d) 

We express all operators in the above equations in terms of the quasiparticle 
operators (A|, Aj,iVj) defined with respect to huiq) and \4>{q)) as follows: 

hM{q) = Y,EiN„ (4-22) 

i 

= {m\ m)) + Y.^a}^^^4 + A) + j;<j(i)Ar„ (4-23) 

i i 

^i"^ = E^i:s^«(A^-^^) (4-24) 



with 



^iti (0 = ^K' - vf), F^^^ii) = 2diaiUiVi, (4-25a) 



^A,/W = -2' ^A:2n0 = 0, (4-25b) 

^B^i (0 = -UiVi, F^j^lii) = dMul - vl)- (4-25c) 



Note that all matrix elements are real so that {4>{q)\ f} "* \4>{q)) = 0. For later 
convenience, we define one-body operators, 

R^s"^ ^ [4^ > ikq) - Kq)N)A] = 2 - a,), (4-26) 

i 

where F^^ is the last terms of Eq. ()4-23|) and 

= {2uMei - xdiOi - \{q)) - A{uj - v^,)] f'^^] {i) . (4-27) 
The infinitesimal generators Q{q) and P{q) can be witten as 

Q{q) = Y.Qi{q){A\ + Ai), (4-28a) 

i 

P(q)=iY,P^mA\-A^)■ (4-28b) 
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Equations ()4-19b|) and (|4-19c|) are then reduced to linear equations for the matrix 
elements Qi{q) and Pi{q) of the infinitesimal generators Q{q) and P{q). They are 
easily solved to give the expression 



Qiiq) 



2Ei 



+ 



s 



(4-29a) 



2E.. 



+ 



s 



i2E,)^ - u;2 



where A'^,- = 2n,t>,- and 



(4-29b) 



(4-30) 



fV^) _ fV-^> _|_ 
J PR,s — JP,S "T JR,s ■ 

Note that a;^, representing the square of the frequency u){q) = ^ B{q)C{q) of the 
local harmonic mode, is not necessarily positive. The values of B(q) and C{q) 
depends on the scale of the collective coordinate q, while u){q) does not. In this 
sense, there remains an ambiguity for determining q. We thus require B(q) = 1 
everywhere on the collective space to uniquely determine q. 

The quantities Jq f^] and are easily calculated to be 

i 
i 

f^^I = 2^sY.Ri^^i^Q^(i)- 

i 

Inserting Eqs. (|4-29ap and H4-29b|) for Qi{q) and Pi{q) into the above expressions, 

we obtain linear homogeneous equations for unknown quantities /p^^ s ' /q i ™^ • 
Similarly, the condition of orthogonality to the number operator, 

{4>iq)\[N,Piq)Mq))=2iJ2NiPiiq)=0, (4-32) 



(4-31a) 
(4-31b) 
(4-31c) 



gives another equation for /p^^^, /g } and /at. Since /g 2 = Oj they are summarized 
in a 4 X 4 matrix form as follows: 



( ' ^ 

A+) 

J PR,2 
\ fN I 



0, 



(4-33) 



where 



Sn=2GS^'-\F^-},F^^})-l, 
5i2 = 2G5(2)(Fi-),Fi+)), 



(4-34a) 
(4-34b) 
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= 2G5(2)(Fi-\F(+)), 


(4-34c) 


Sii 




(4-34d) 


S21 




(4-34e) 


S22 




(4-34f) 


S23 


= 2G{5W(FW,Fl+)) + 5(^)(i?W,Fi+))}, 


(4-34g) 


S24 


= 2G{S'(i) , N) + 5(2) (4+) , AT)} , 


(4-34h) 


S31 




(4-34i) 


S32 




(4-34j) 


S33 


= ^x{S^'HFilF^^^) + Si^HR^^\Fi^^)} - 1, 


(4- 34k) 


S34 


= 4x{5(i) , TV) + 5(2) (4+) , TV)}, 


(4-341) 


S41 


= u;^5(^)(Ar,Fi-)), 


(4-34in) 


S43. 


= 5W(iV,Fi+)), 


(4-34n) 


Sa3 


= 5W(iV,FW), 


(4-34o) 


Saa 


= 5W(iV, iV). 


(4-34p) 



Here, the functions S^^\X,Y) and S''^\X,Y) are defined by 

i 

S('\X,Y) = ^ ^^2__x,y„ (4.35b) 

with Xi and denoting ones among the quantities F^^ {i), Ri'^\i) and Ni. The 
frequency a;(g) is determined by finding the lowest solution of the dispersion equation 

det{Skk'iuj)} = 0, (4-36) 

i.e., the solution of which lo'^ is minimum. Normalizations of /p^\, /p/j2' /q 1 
/jv are fixed by the canonical variable condition 

(<^(3)|[Q(g),F(g)]|<^(g)) =2i^Qi(g)P,(g) =i. (4-37) 

i 

Note that represents the curvature of the collective potential, 

for the choice of coordinate scale such that B{q) = 1. 

§5. Numerical analysis 



In this section, we make numerical analysis of the oblate-prolate shape coexis- 
tence and large-amplitude collective motions in the multi-0(4) model. 
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5.1. Procedure of calculation 

We first solve the HB equations and find HB equilibrium points which corre- 
spond to extrema of the collective potential V{q) defined by Eq. (|2-12|) . At these 
points, the HB equation in the moving frame, Eq. ()2-7() . and the local harmonic 
equations, Eqs. (|2-9a() and (|2-9b|) . coincide with the ordinary HB equation and the 
quasiparticle RPA equations, respectively. Let |(/)((?o)) be a HB solution, which is 
assumed to be on the collective subspace at a particular value of the collective coor- 
dinate, q = qo. Solving the quasiparticle RPA equation with respect to \(j){qo)), we 
find a collective normal mode, which determines the infinitesimal generators Q{qo) 
and P{qo). In the present analysis, we choose the normal mode with the lowest 
frequency, which is most collective with respect to to the "quadrupole" operator D. 
We then generate the state |(/'((?o + ^q)) with an infinitesimal shift of the collective 
coordinate as 

|^(go + <5g))=e-^^^(^")|c/>(go)). (5-1) 

Next, we solve the local harmonic equations at q = qo + 6q and determine Qiqo + ^Q) 
and P{qo + 6q), and proceed to q = qo + 26q. Repeating this procedure, one can 
construct a collective subspace. Thanks to the invariance, 

4 -Ai = a\ - Au (5-2) 

we can rewrite the state vectors as 

m + Sq)) = exp{dg ^P,(g)(Al - A,)} |</>(g)) 

i 

= exp{5g^P,(g)(4 - A,)} \cPiq)) . (5-3) 

i 

Combining the above expression with Eq. (|4-1|) . we obtain the following simple 
relations between the Ui,Vi coefficients at q + 6q and those at q: 

f Ui{q + 6q) \ ^ f cos{Pi{q)Sq) -sm{Pi{q)6q) \ / Ui{q) \ , , 

\vi{q + 6q)J \ sm{Pi{q)6q) cos{Pi{q)6q) ) \ Vi{q) ) ' ^> 

In the present calculation, we always start from the spherical equilibrium point 
{q = 0). This point is an unstable extremum (a saddle point) on the collective 
subspace when we choose parameters producing deformed HB minima (see below). 
Note that the quasiparticle RPA equations given in subsection 4.1 are still valid 
whereas the frequency uj of the eigenmode is pure imaginary in this case. 

If the collective subspace is exactly decoupled from the non-collective subspace, 
the state vectors \4'{q)) obtained in this way should satisfy the HB equation in 
the moving frame, Eq. ()2-7|) . simultaneously. In general, however, the decoupling 
conditions may not be exactly satisfied and we have to resort to some kind of iterative 
procedure to construct a collective subspace which complies with Eq. l|2-7() . One way 
of evaluating the quality of decoupling and accuracy of the numerical calculation is 
to examine validity of Eq. ()2-7|) on the collective space. For the multi-0(4) model 
under consideration, as we shall see in Fig. 2, this condition is found to be well 
satisfied with the use of a step size 5q = 0.005 in Eq. H5-4|l . Of course, this does 
not necessarily mean that such a simple algorithm will always work also for more 
realistic cases, and we will need more investigations concerning numerical techniques 
of solving the set of equations (|^ . ^F^, and (|2-9b|l . 
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The collective Hamiltonian thus obtained, 'H{q,p) = ^B{q)p^ + V(q), is then 
quantized, and the collective Schrodinger equation is solved to get eigenvalues and 
transition probabiHties. As wc set a scale of the collective coordinate q such that 
B(q) = 1, there is no ambiguity of ordering in the canonical quantization procedure, 
following the Pauli's quantization rule. It is easily confirmed that the collective repre- 
sentation of the "quadrupole" operator, defined by X'(g) = {4>{q,p) \ D \ (l){q,p)), does 
not depend on p, and transition matrix elements can be evaluated by / ■0n(Q')*^(g)'0n' {Q)dq, 
where ipniQ) denotes the collective wave function of the n-th eigenstate. 

5.2. Parameters 

In the following, wc consider a case consisting of three shells with the spherical 
single-particle energies, e^^ = 0.0, e^^ = 1.0, e^^ = 3.5, the pair degeneracies, 
f2j^ = 14, i7j2 = 10, f2jg = 4, the reduced quadrupole moments, dj-^ = 2, dj^ = 
dj^ = 1, and we distribute 14 pairs {Nq = 28) in this shell-model space (wc do not 
distinguish protons and neutrons). We compare numerical results obtained with 
different values of the pairing interaction strength, G = 0.14, 0.16, and 0.20, for a 
fixed "quadrupole" interaction strength % = 0.04. These numerical examples are 
presented merely as representatives of similar results obtained with other sets of 
parameters. Since properties of the single-shell 0(4) model are determined by the 
ratio G/x, we obtain similar results varying x instead of G. The only reason why 
we here vary G fixing % is the convenience of visualizing the change of the barrier 
height (between the prolate and oblate local minima) while keeping the magnitude 
of the equilibrium deformation Dq almost constant. 

Figure 1 illustrates the single-particle energy diagram as a function of D. Details 
of the level crossings are explained in the caption to this figure. It is certainly 
possible to discuss the level crossing dynamics and simulate the shape coexistence 
phenomena by means of the multi-0(4) model consisting of only two shells. However, 
consideration of the three shells is found to be more appropriate in order to make 
properties in the large deformation region more realistic. In realistic situations, a 
number of successive level crossings take place with increasing deformation. 

5.3. Collective potentials 

We present in Fig. 2 the collective potentials calculated for three values of the 
pairing strength G. For G = 0.14, it exhibits a double well structure, while it looks 
like an anharmonic-oscillator potential for G = 0.2. The case of G = 0.16 shows 
a transitional character between these two situations. Note that the collective po- 
tentials arc plotted here as functions of deformation D. The relation between D 
and the collective coordinate q will be explicitly shown later. The mechanism of 
appearance and disappearance of the double well structure is determined by compe- 
tition between the pairing and "quadrupole" interactions: The pairing force favors 
an equal population of particles over all magnetic substatcs in each j-shell, whereas 
the "quadrupole" force favors an occupation of down-sloping levels creating level 
crossings and asymmetry in the magnetic substate population. In the three j-shell 
case under consideration, the energy of the degenerate local minima is lower than 
that of the spherical point for small values of G. With increasing pairing correla- 
tions, however, the energy of the spherical point decreases, so that the barrier hight 
decreases and eventually the barrier itself diminishes. 

In Fig. 2, quantum eigenstates obtained by solving the collective Shrodinger 
equations are indicated. We can clearly see a parity-doublet-like pattern associated 
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Fig. 1. Deformed single-particle energies, e; = e^—xdiOiD, plotted as functions of the deformation 
parameter D for a set of parameters: ej^ — 0.0, e^^ = 1-0, = 3.5, flj-^ — 14, Qj^ = 
10, i7j3 = 4, and dj-^ = 2, dj^ = dj.^ = 1. Numbers enclosed by circles indicate the numbers 
of pairs occupying individual levels in the lowest energy configuration for a given value of D 
when pairing correlations are absent. Near the spherical point [D — 0), the lowest shell is 
fully occupied by 14 pairs (28 nucleons). The level diagram is symmetric with respect to the 
sign-change of D (prolate-oblate symmetry). The down(up)-sloping levels (with increasing D) 
are named "dl,d2, and d3 (ul,u2, and u3)". In the prolate side {D > 0), the down-sloping 
level originating from the second shell (d2) crosses the up-sloping level (ul) at D = 8.3. After 
this first crossing, the lowest-energy configuration involves 5 pairs occupying the second down- 
sloping level (d2) and 2 pairs remaining in the first up-sloping level (ul). At D = 29.2, the 
down-sloping level originating from the third shell (d3) crosses the second up-sloping level 
(u2). In the lowest-energy configuration after this crossing, all the down-sloping levels are 
fully occupied while all the up-sloping levels are unoccupied. This configuration possesses the 
maximum value of the deformation parameter, Dmax ~ 42. The same pattern of level crossings 
is valid also for the oblate side {D < 0). Note that this figure is drawn by regarding D as a 
free parameter, although D in Eq. H4-5II means an expectation value of the operator D. In this 
sense this figure should be regarded as a kind of Nilsson diagram. 

with tunneling through the central barrier in the collective potential curve. Prop- 
erties of these excitation spectra will be discussed in the succeeding subsection. In 
this figure, we also show the potential energy functions calculated by means of the 
conventional constrained HB (CHB) procedure: 

K:hb(I?) = (<A?™(Z))|i?|<?^?™(Z))) = 2j2e^vf - ^ - '-xD', (5-5) 

i 

where the state vectors \(f)Q^^{D)) are determined by the constrained variational 
principle 

5(</><r™(Z))|i? -XN- f,D\cl>^^^m = 0, (5-6) 

with /i denoting a Lagrange multiplier. As we see, the collective potentials obtained 
with the ASCC and CHB methods are practically indistinguishable. In principle, 
when both the collective path obtained in the ASCC method and that in the CHB 
method go through the same HB local minima, the collective potential energies 
at these points should coincide with each other. One may notice, however, very 
small differences between their values at the HB local minima with D ^ for the 
cases of G = 0.14 and 0.16. These discrepancies are due to violation of Eq. 1)2 -711 
which accumulates in the numerical calculation starting from D = 0, and indicate 
the amount of error associated with the computational algorithms adopted here, as 
mentioned in subsection 6.1. 
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Fig. 2. Collective potentials calculated with the ASCC method for G = 0.14 (top), G = 0.16 
(middle), and G = 0.20 (bottom) arc plotted by solid linos as functions of deformation D. 
Other parameters used in the calculation are listed in the caption to Fig. 1. For comparison, 
collective potentials of the CHB method are indicated by dotted lines. Horizontal solid lines 
indicate eigen-energies of the quantized collective Hamiltonian. 

In Fig. 3, the pairing gaps A are shown as functions of deformation D. They 
monotonically decrease as D increases, and vanish at the maximiim deformation, 
which is Dmax = 42 for the parameters adopted in these calculations. In Fig. 4, 
occupation probabilities vf are displayed. Note that all down-sloping (up-sloping) 
levels are fully occupied (unoccupied) at i^max- Again, the results obtained with the 
ASCC and CHB methods are practically indistinguishable in these figures. Thus, 
in the present multi-0(4) model, the static mean-field properties, like the potential 
energy curves and the pairing gaps, obtained with the ASCC method, are very 
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similar to those obtained with the conventional CHB method. This is due to the 
simplicity of the multi-0(4) model adopted here, and the collective subspace in 
general may differ from that of the CHB method I27|l . 
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Fig. 3. Pairing gaps A calculated with the ASCC method for G = 0.14 (top), G = 0.16 (middle), 
and G — 0.20 (bottom) are plotted by solid lines as functions of deformation D. Other param- 
eters used are the same as in Fig. 2. For comparison, A calculated with the CHB method and 
with the BCS approximation are also shown by dotted and dashed lines, respectively. The solid 
and dotted lines appear similar to each other. The equilibrium deformations are indicated by 
arrows. 



In Figs. 3 and 4, we have also shown the BCS pairing gaps and occupation prob- 
abilities calculated in the BCS approximation. In this approximation, single particle 



energies given by 



xdiGiD are used to solve the BCS equations, where the 



deformation D is treated as an ad hoc parameter, neglecting the HB selfconsistency 
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Fig. 4. Occupation probabilities v\ calculated with the ASCC method for G = 0.14 (top), G — 0.16 
(middle), and G = 0.20 (bottom) are plotted by solid lines as functions of the deformation D. 
Other parameters used are the same as in Fig. 2. For comparison, vf calculated with the 
CHB method and with the BCS approximation are also shown by dotted and dashed lines, 
respectively. The solid and dotted lines appear similar to each other. Note that, at the large 
deformation limit Dmax, all down-sloping (up-sloping) levels are fully occupied (unoccupied) in 
the cases of the ASCC and CHB methods, whereas it is not the case in the BCS approximation. 



H4-5I) in accord with phenomenological singe-particle potential models. The pairing 
gap A in the BCS approximation differs from those in the CHB and ASCC methods. 
Note that the BCS gaps do not vanish in the limit i^max in contrast to the HB self- 
consistent pairing gaps. In the following we shall focus our attention on dynamical 
properties of tunneling motions between the oblate and prolate minima. 
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5.4. Excitation spectra and collective mass 

Figure 5 shows the excitation spectra and transition matrix elements for the 
ASCC method and for the exact diagonahzation in a basis set, Eq. (jSIHl)- Wave 
functions of the low- lying states are displayed in Figs. 6 and 7. The collectivity of 
these states are apparent from the enhancement of the transition matrix elements 
in comparison with the single-particle ones dj{=l or 2). We see that the ASCC 
method well reproduces the major characteristics of the exact spectra and transition 
properties. In view of the huge degrees of freedom involved in the model under 
consideration (the dimension of this shell model space is 1894), it is remarkable that 
the low-lying states can be described very well in terms of only the single collective 
coordinate q. In particular, we note that the emergence of the "parity splitting" 
pattern for weaker strength G of the pairing interaction is well reproduced in the 
calculation. This implies that the ASCC method succeeds in describing the large 
amplitude tunneling motions through the barrier between the oblate and prolate 
local minima. 

As is well known, the collective mass parameter represents inertia against change 
of the mean field. It is determined locally and changes as D varies. Since we have 
set the scale of the collective coordinate q such that the collective mass parameter 
M{q) = B^^{q) with respect to q is unity, the collective kinetic energy can be written 
either by g or D as 

Thus, we obtain the explicit expression for M{D): 

M{D) ^ {^f = (±Y.2d,a,vKq))-' = {i^dia.umPMr' ■ (5-8) 

i i 

The collective mass M{D) evaluated in this way is shown in Fig. 8 as a function of D. 
One immediately notice that the ASCC mass M[D) diverges in the limit D I?max- 
The reason is clearly understood by examining the relationship between D and the 
collective coordinate q. It is shown in Fig. 9. We see that the value of D saturates as 
it approach its limiting value -Dmax- Thus, the derivative dq/dD, which corresponds 
to the collective mass M{D) according to Eq. (|5-8|) . diverges. Needless to say, this 
divergence is caused by the existence of the maximum value of D which is an artifact 
of the model: If we increase the number of shells explicitly taken into account, this 
limit is removed by successive level crosings with increasing D. 

It is interesting to compare the ASCC mass M{D) with the cranking mass 



En{D)-EoiD) 

\2uMxdiai + §^) + {uj-vf)^\'' 
{2Eif 



(5-9) 



where 4>o^^{D) and represent the ground and excited states, and Eq{D) 

and En{D) the energies of them obtained in the BCS approximation. Namely, the 
coefficients of the Bogoliubov transformations, Ui, Vi, the pairing gap A, the chemical 
potential A, and the quasiparticle energies Ei are evaluated with use of single-particle 
energies defined by = e^ — xdiOiD. It is important to note that the deformation D 
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Fig. 5. Comparison of excitation spectra calculated with the ASCC method and and with the 
exact diagonalization for G — 0.14 (top), G = 0.16 (middle), and G = 0.20 (bottom). For 
reference sake, excitation energies of the lowest RPA modes in the HFB local minima are also 
indicated. Other parameters used are the same as in Fig. 2. Numbers adjacent to vertical 
lines indicate transition matrix elements for the "quadrupole" operator D. In the top panel, 
transition matrix elements between "parity doublets" are indicated with arrows. 



is treated in the BCS approximation as an ad hoc potential parameter disregarding 
the selfconsistency condition (|4-5j) . Figure 8 shows that the cranking mass Mcr{D) 
is significantly different from the ASCC mass M{D) in the whole interval of the 
deformation D, including the spherical point D = and the deformed equilibrium 
points. The difference between the ASCC mass and the cranking mass can be 
understood in terms of the HB selfconsistency. 

At the spherical point (-D = 0), we can make the comparison between the ASCC 
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Fig. 6. Collective wave functions -tpiD) for the ground and first excited states in Fig. 5, obtained 
with the ASCC method for G = 0.14 (top), G = 0.16 (middle), and G = 0.20 (bottom). The 
equilibrium deformations are indicated by arrows. The ipiD) are defined by i){D) = i'{g)\-§^\'^ , 
where tp{q) are wave functions for the collective coordinate q. Thus, they are normalized as 
/ \i^{q)\'dq = f \i^{q)\'\§§\dD = f mD)\^dD = 1. 



mass and the cranking mass in an explicit way. There, all terms linear with respect 
to (Tj in the local harmonic equations vanish after summing over all levels i so that 
the pairing and "quadrupole" normal modes are exactly decoupled. Thus, we get a 
simple expression of the ASCC mass 

M(D = 0) = 2x' T = Ax'A' T ^ (5-10) 
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Fig. 7. The same as Fig. 6 but for the second and third excited states. For G — 0.2, only 
the wave function for the second excited state is drawn. 



On the other hand, the expression of the cranking mass Mct{D) reduces to 

i 



(5-11) 



since the derivatives, dA/dD and dX/dD vanish at D = 0. We see that the above 
expression for Mcr{D = 0) can be related to M{D = 0) if we set (J = in Eq. (f^THjl . 
In reality, the frequency uj is imaginary and uj'^ is negative when the barrier exists 
{G = 0.14 and 0.16). Accordingly, M{D = 0) < M^riD = 0) in these cases. On the 
other hand, w is real and is positive when the spherical point is stable {G = 0.20), 
so that M[D = 0) > Mcr{D = 0) in this case. In this way, the difference between 
the ASCC mass and the cranking mass in the barrier region (near D k. {)) noticeable 
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Fig. 8. Collective mass M{D) calculated with the ASCC method for G = 0.14 (top), G = 0.16 
(middle), and G — 0.20 (bottom) are plotted by solid lines as functions of deformation D. Other 
parameters used are the same as in Fig. 2. For comparison, the cranking mass are indicated by 
dotted lines. The equilibrium deformations are indicated by arrows. 

in Fig. 8 can be understood in terms of the finite frequency effect of the local 
harmonic mode, which is connected with the curvature of the collective potential 
by Eq. ()4-c{8() . This finite frequency effect decreases (increases) the collective mass 
when the barrier emerges (diminishes). We would like to emphasize that the finite 
frequency effect in the local harmonic equations represents selfconsistent dynamics 
of the time-dependent mean-field. 

We note that, in contrast to the ASCC mass M{D), the cranking mass Mcr{D) 
does not diverge at Dmax- This is because, as already emphasized above, the de- 
formation D is treated as an ad hoc parameter in the single-particle potential, so 
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Fig. 9. Relations between the deformation D and the coUective coordinate q in the ASCC 
method. The soHd, dashed and dotted hnes show the results of calculation for G = 
0.14,0.16, and 0.20, respectively. 

that the existence of the hmit -Dmax in the selfconsistent deformation defined by 
Eq. H4-5() is disregarded there. In the multi-0(4) model under consideration, the 
selfconsistency between the deformation parameter specifying the single-particle po- 
tential and the density deformation evaluated in terms of the wave function becomes 
extremely important near -Dmax- 

The comparison between the ASCC and cranking masses in Fig. 8 thus indicates 
the importance of the HB selfconsistency in evaluating the collective mass, although 
it is somewhat exaggerated because of an artificial aspect of the present parameter 
setting of the multi-0(4) model. 

There are various microscopic approaches to derive the collective mass (also 
called inertial functions) .^''^''^ Certainly, it is important and interesting to make 
a detailed comparison of different approaches for the collective mass in the multi- 
0(4) model, and clarify the role of the HB selfconsistency in more detail. Such a 
more systematic and comparative study is beyond the scope of this paper, but will 
be pursued in a separate paper. 

§6. Conclusions 

The ASCC method was applied to the exactly solvable multi-0(4) model which 
simulates nuclear shape coexistence phenomena. The collective mass and dynamics 
of large amplitude collective motions in this model system were analysed, and it was 
shown that the method can well describe the tunneling motions through the barrier 
between the prolate and oblate local minima in the collective potential. The result 
of numerical analysis strongly encourages applications of this approach to realistic 
cases. We plan to investigate the oblate-prolate shape coexistence phenomena in 
^^Se with use of the P-l-Q interactions. 
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